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Abstract. A new approach to the finite-gap property for the Heun equation is 
constructed. The relationship between the finite-dimensional invariant space and 
the spectral curve is clarified. The monodromies are calculated and are expressed as 
hyperelliptic integrals. Applications to the spectral problem for the BC\ Inozemtsev 
model are obtained. 



The BCi Inozemtsev model is a one-particle quantum mechanics model whose 
Hamiltonian is given as 



where p(x) is the Weierstrass p-function with periods (l,r), ou — 0, U\ — 1/2, 10% = 
— (r + 1)/2 and oj^ = r/2 are half-periods, and k {i = 0, 1, 2, 3) are coupling constants. 
This model is a one-particle version of the BCn Inozemtsev system [7j, which is known 
to be the universal quantum integrable system with Bn symmetry The BCn 

Calogero-Moser-Sutherland systems are special cases of BCn Inozemtsev systems. 

It is known that there is a correspondence between the BC\ Inozemtsev model 
and the Heun equation. Here the Heun equation is a standard canonical form of the 
Fuchsian equation with four singularities defined in (Jo"7|) (see section 

The potential Yui=o ^(^ + + u i) nas a l so arisen from a different source, the 

theory of elliptic solitons. Around 1987-92 Treibich and Verdier JU| found and showed 
that Y^h=o + + u i) 1S an algebro-geometric finite-gap potential iff li G Z for 

all i G {0, 1, 2, 3}. For this reason the function ^2 i=0 k(k + l)p(^ + ^>i) is called the 
Treibich- Verdier potential. Subsequently several others have produced results on this 
subject, e.g., Gesztesy and Weikard jU |2Ij, Smirnov [TJ] and so forth. (For more 
details see the references in |H |2*T] ITij.) 

In this paper a new approach to the Treibich- Verdier potential is presented, along 
with applications to the spectral problem for the BC\ Inozemtsev model. In particular 
the relationship between some finite-dimensional space of elliptic functions and the 
spectral curve is clarified and a new determinantal formula for a commuting operator 
of odd degree is obtained. Here the existence of a commuting operator of odd degree 
is equivalent to the algebro-geometric finite-gap property. 

The reader may find the author's background in this subject to be helpful in eluci- 
dating the motivations for this research. In ^H] I proposed and developed the Bethe 
ansatz method for the Heun equation and applied it to the analysis of the eigenstates 
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and the eigenvalues for the BC\ Inozemtsev model. A few months after the work ^H] 
I noticed that the function ^i=o ^« + l)fp( x + w i) ls an algebro-geometric finite-gap 
potential while still being unaware of the results of Treibich and Verdier. An argu- 
ment based on the Bethe ansatz led me to conjecture that the global monodromies 
are expressible in terms of hyperelliptic integrals and that they would be applicable 
for investigating the singularities of the eigenvalues in r. Eventually the results of 
Treibich and Verdier [T^j, Hochstadt jH], McKean and van Moerbeke ^U] and Smirnov 
[T4"] came to my attention. In particular some formulas for the monodromy are given 
in jHl [TUJ and it seemed natural that they would be applicable for the Treibich- Verdier 
potential. 

We next present a brief sketch of the contents and results of this paper. 

Firstly a method of determining whether the function Y^=oh(h + + ^>i) 

(lo,h,h,h £ Z) is an algebro-geometric finite-gap potential is presented in sections 
E] and El The strategy is as follows: First the product of some two solutions to 
H f{x) = Ef(x) (with H the Hamiltonian, E the eigenvalue) is shown to be doubly 
periodic for all E, a result which was already obtained in ^B] i n order to show the 
completeness of the Bethe ansatz. From the solution to the equation for that product 
a differential operator A of odd degree that commutes with the Hamiltonian H is 
constructed in section |3~T1 whence follows the algebro-geometric finite-gap property. 
The proof is elementary and does not require the advanced techniques of Treibich and 
Verdier that are based on the theory of Jacobi varieties and the notion of a tangential 
cover. 

In section 12.21 finite-dimensional spaces of elliptic functions are introduced. Note 
that these finite-dimensional spaces have already appeared in |1& and they are related 
to the concept of quasi-exact solvability [201 E] These spaces play an important role 
in the presentation and investigation of the spectral curve, i.e. the curve determined 
by the equality between the Hamiltonian H and the commuting operator A. As an 
application we obtain a determinantal formula for the commuting operator in section 
13.21 In section 13.31 formulas for the monodromy of the solutions to the Inozemtsev 
model are obtained. These formulas are expressed in the form of hyperelliptic in- 
tegrals. In section 13.41 examples are used to illustrate the results and the formulas 
obtained in sections 121 and El 

From the viewpoint of physics it is often desirable to obtain the square-integrable 
eigenfunctions of the Hamiltonian and the corresponding eigenvalues. In sections 01 
and E3 the eigenfunctions of H without this physical boundary condition are con- 
structed, whereas in section HI the physical boundary condition is included. 

It is shown in [TTj that an eigenvalue satisfying the physical boundary condition is 
holomorphic in r if —1 < p = exp(ny/— It) < 1. In this paper a sufficient condition 
is provided for the holomorphy of the physical eigenvalues in r G C. This condition 
involves the convergenec radius of the eigenvalue as a power series in p, which can be 
calculated by the method of perturbation [T7| . 

In section |4~T1 the strategy is outlined. Section l4~2l provides the appropriate setting 
for consideration of the spectral problem in the case lo,h G Z>i. In particular the 
Hilbert space H and its dense subspace D are introduced. In section 1431 we obtain 
a sufficient condition for the holomorphy of the physical eigenvalue in r G C for this 
case. Here the formulas for the monodromy from section 13.31 are used. An example 
is also presented to illustrate the results. In section l4~4l the case l G Z>i and l\ = 
is discussed and in section 14.51 the case Iq = l\ = is treated. 
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In section |3] the relationship between the Heun equation and the Inozemtzev model 
is reviewed and a translation between the two formalisms is detailed. Formulas for 
the monodromy along the cycle around two regular singular points are derived. They 
are expressed in use of hyperelliptic integrals. 

In the appendix definitions and formulas for the elliptic functions are provided. 

Throughout this paper it is assumed that U G Z {i = 0, 1, 2, 3) and (Zo, h, h, h) 7^ 
(0,0,0,0). 



2. Results from [16J and related topics 

2.1. Product of the solutions. 

Fix the eigenvalue E of the Hamiltonian H and consider the second-order 
differential equation 

(2) (H - E)f(x) = + + l)p(* + Wi) " E^j f(x) = 0. 

Let h(x) be the product of any pair of the solutions to (J2J). Then the function h(x) 
satisfies the following third-order differential equation: 

(3) 

^ - 4 (j2 + l M x + w *) - ^ ^ - 2 [Y. k{li + + J h ^ = °- 

It is known that equation (J3J) has a nonzero doubly periodic solution for all E if 
li G Z>o (i = 0, 1, 2, 3). That was first observed by Darboux pQ. 

Proposition 2.1. l?Ti Proposition 3.5] If lo,h,h,h G Z> , i/ien equation /ias a 
nonzero doubly periodic solution S(x,E), which has the expansion 

3 k-1 

(4) Z(x,E) = ^E) + J2Y, b ?( E ^ x + UJ ^ J > 

i=0 j=0 

where the coefficients cq(E) and b^\E) are polynomials in E, they do not have com- 
mon divisors and the polynomial cq(E) is monic. We set g = deg E cq(E) . Then the 
coefficients satisfy deg E by (E) < g for all i and j . 

We can derive the integral formula for the solution A(x, E) to the equation (J2J) in 
terms of the doubly periodic function E), which is obtaind in ^H]- Set 

(5) 

Q(E) = ~(x, Ef (e-Y, U{U + l)p(x + ) + \e(x, E)^§^ - ~ 

It is shown in that Q(E) is independent of x. Thus Q(E) is a monic polynomial 
in E of degree 2g + 1, which follows from the expression for i?) given by (j^J). The 
following proposition is obtained in |16j : 

Proposition 2.2. [TBJ Proposition 3.7] Let S(x,E) be the doubly periodic function 
defined in Proposition \2. 1\ and Q{E) be the monic polynomial defined in ^j). Then 
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the function 

'-Q{E)dx 



(6) A(x,E) = y^(x~E) e x P J ^_ 

is a solution to the differential equation 



2.2. The invariant subspace of elliptic functions. 

It is shown in ^Hj that there are finite-dimensional spaces of elliptic functions 
on which the action of the Hamiltonian JJJ) is well defined. In this subsection, we 
introduce such spaces and discuss the relationship with the polynomial Q(E) of the 
previous subsection. 

Let T be the space spanned by meromorphic doubly periodic functions up to signs, 
namely 

ei,e 3 =±l 

(8) T t ^ 3 = {f{x)\ meromorphic \ f(x + 2ui) = e\f{x), f(x + 2u 3 ) = e 3 f(x)}, 

where (2a>i, 2u 3 ) are basic periods of elliptic functions. 

Let U G Z>o (i = 0,1,2,3) and let ki be the rearrangement of U such that ko > 
k\ > k 2 > k 3 (> 0). The dimension of the finite-dimensional invariant subspaces of 
the space T and the construction of the spaces are given in ^Hi Theorem 3.1]. 

Proposition 2.3. [TBI Theorem 3.1] 

(i) If the number k + ki + k 2 + k 3 {= l + l x + l 2 + l 3 ) is even and k + k 3 > k x + k 2 , 
then the dimension of the maximal finite- dimensional invariant subspace in JF with 
respect to the action of the Hamiltonian H (see ^) is 2k + 1. 

(ii) If the number k + k\ + k 2 + k 3 is even and k + k 3 < k\ + k 2 , then the dimension 
is k Q + ki + k 2 -k 3 + 1. 

(Hi) If the number k + k\ + k 2 + k 3 is odd and k > k± + k 2 + k% + 1, then the 
dimension is 2ko + 1. 

(iv) If the number ko + k\ + k 2 + k 3 is odd and ko < ki + k 2 + k 3 + 1 , then the 
dimension is k + k\ + k 2 + k 3 + 2. 

Let V be the finite-dimensional invariant space described in Proposition 12.31 and 
P{E) be the monic characteristic polynomial of the Hamiltonian H (JJJ on the space 
V, i.e. P(E) = US V ( E ~ E%) {{Ei} are eigenvalues of H on V). 

Here we present the invariant space V more explicitly. We assume G Z (i = 
0,1,2,3) and -EtoA/2 G Z> . Let 

AAA ^ e ^ e ^ l)-dimensional vector 
space spanned by < pi(a;) /3l p 2 (x) /32 p3(x)^ 3 p(x) ?1 i , where d = — Ei=o A/2 and 

pi(x) (z = 1,2,3) are co-p functions (see appendix). We assume on G {—k,h + 1} 
(z = 0, 1, 2, 3) and Eto on/ 2 G Z. Then we set 

{Vao,ai,a2,a3i Ej=0 ^i/2 ^ Z<q, 

{0}, Eto«i/ 2 = 1 ! 

^l-a ,l-ai,l-a2,l-a 3 ) 
Eto«i/ 2 eZ> 2 . 

If /o + ^i + ^2 + ^3 is even, then the space V is decomposed as V = £/_j 0) _j 1) _fe ) _j 3 © 
^-io-Ji,«2+i,i3+i©^-«QA+i-«2,i3+i©^-io,ii+i,i2+i-«3- If + + + ^3 is odd, then we 
have V" = U-i -i lt -i 2t i 3+ i © £^-z ,-2i,z 2 +i,-z 3 © £/-i ,/i+i,-k,-*3 © l^+i,-^ ,-i 2 ,-j 3 . 
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For the equations P(E) = and Q(E) = we obtained the following proposition 
in HE|. 

Proposition 2.4. (c.f. PHI Theorem 3.8]) The set of zeros of Q(E) coincides with 
the set of zeros of P(E). 

We propose the following conjecture. 

Conjecture 1. P(E) = Q(E). 

By applying Proposition 12.41 we immediately have the following proposition. 

Proposition 2.5. (i) If there are no multiple roots for each of the equations P(E) = 
and Q(E) = for generic r, then Conjecture^ is true. 

(ii) If there are no multiple roots for the equation Q(E) = for generic r and 
deg E P(E) = deg E Q(E), then Conjecture^ is true. 

The following proposition follows upon combining Theorem 3.2 and Proposition 
3.9 in ^1 with Proposition 12.51 (ii) of this paper. 

Proposition 2.6. // two of the U (i — 0, 1, 2, 3) are zero, then Conjecture^ is true. 

If Z + 1\ + Z 2 + h is small, then we can confirm Conjecture directly through direct 
calculation. 

Proposition 2.7. (c.f. [TBJ Proposition 3.10],) Iflo + h + h + h < 8, then Conjecture 
Q] is true. 

3. Finite-gap property 

3.1. The construction of the commuting operator. 

Let us consider the Schrodinger operator —d 2 /dx 2 + q(x) with the conditions that 
q(x) is meromorphic on R (may have poles) and q(x + 1) = q(x). These conditions 
on q(x) are assumed in the definitions below. 

Definition 1. Assume that q(x) G M for x G 1R and there are no poles for q(x) on IR. 
We set H = —d 2 /dx 2 + q(x). Let o-b(H) be the set such that 

E G a b (H) <^> All solutions to (H - E)f(x) = are bounded 01116R 

and o~b{H) be the closure of the set o~b(H). If the set R \ cr^H) can be written as 

(10) R \ a b (H) = (-00, Eq) U(E 1 ,E 2 )U---U (£ 2g _i, E 2g ) 

with Eq < Ei < ■ • ■ < E 2g , then q(x) is called the finite-gap potential. 

Definition 2. If there exists an odd-order differential operator A = (d/dx) 29+1 + 
Y^=q bj(x) (d/dx) 29 ~ 1 ~ : ' such that [A, —d 2 /dx 2 + q(x)] = 0, then q(x) is called the 
algebro-geometric finite-gap potential. 

Remark. It is known that the equation [A, —d 2 /dx 2 + q(x)] = is equivalent to the 
function q(x) being a solution to some stationary higher-order KdV equation. 

It is known that, if q(x) is real-holomorphic on R and q(x + 1) = q(x), then q(x) is 
a finite-gap potential if and only if q(x) is an algebro-geometric finite-gap potential 
(for details see [UJiaini)- 
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In this section it is shown that Yli=o^i(^ + l)fp( x + u i) is an algebro-geometric 
finite-gap potential if U e Z (z = 0, 1,2,3). By the substitution /j «-> — Zj — 1, it is 
sufficient to show it for the case U G Z> (i = 0, 1, 2, 3). 

Write 

9 

(11) E(x,E) = J2^(x)E\ 

i=0 

From Proposition 12. II we have a (x) = 1. 

Theorem 3.1. Set u(x) = X^=o^(^ + l)fp( x + w i) anc ^ define the (2g + l)st-order 
differential operator A by 

< 12 ) * = f {<*<*>£ ~ \ (^W) } (-£ + «(»))' ' > 

where the dj(x) are defined in Ml)) . Then the operator A commutes with the Hamilton- 
ian H = —d 2 /dx 2 + u(x). In other words the function u(x) = Yl^=o h{h + l)p( x + u i) 
is an algebro-geometric finite-gap potential. 

Proof. Since the function E(x,E) in (fTTj) satisfies the differential equation Q, we 
obtain the following relations by equating the coefficients of E 9 ~i\ 

(13) a"'(x) — 4u(x)a'j(x) — 2u'(x)aj(x) + 4a^ +1 (x) = 0. 

On the other hand 

9 r 



(14) [A,H] = J2 



3=0 



dx 



H g-j 



:a!"{x) -u'(x)a j (x) ) H 9 ~ j 



3=0 v 

(~mx)(-H + u(x)) + ~af(x) - u , (x)a j (x)^j H 9 ^ 



3=0 

9 



\ ( 4a i+i( x ) - M^Vjix) + Or) - 2u'(x)a j (x)) H 9 ~ j 



2 

j=0 



Proposition 3.2. Let Q(E) be the polynomial defined in (0). Then 

(15) A 2 + Q(H) = 0. 

Proof. Let f(x) be a solution to the second-order differential equation (H — E)f(x) 
0. From the definitions of A in (J12)) and (Jll)) it follows that 

(16) Af(x) = £ { aj (x)E'-'± - \ (jL aj (x)Er-^ } f(x) 

= ~(x,E)f(x)-^~'(x,E)f(x). 



□ 
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Here E'(x,E) := J^E(x,E). From the commutativity of A and H we have (H — 
E)Af(x) = 0. Hence 

(17) A(Af(x)) = E(x, E)± (Af(x)) - ±E'(x, E)Af(x) 

= E(x,E) 2 f(x) - ~E"(x,E)E(x,E)f(x) + \e\x, Eff{x). 

= (e(x, Ef{u{x) -E)- l -E"{x, E)E(x, E) + ±E'(x, E) 2 ^j f(x) 
= -Q(E)f(x) = -Q{H)f{x), 

where relation (j^J) was used. The operator A 2 + Q(H) is a differential operator with 
degree at most 2g + 1. Since the equality (A 2 + Q(H))f(x) = is valid for all 
eigenfunctions of H, it follows that A 2 + Q(H) =0. □ 

Remark. Theorem Kill and Proposition 13.21 are established in for a smooth 
finite-gap potential. 

Let r : y 2 = —Q(x) be the spectral curve of the operators A and H. This curve 
plays an important role in the formulas © and ()15j) . By Proposition 12 .41 edges of the 
hyperelliptic curve T are eigenvalues on the invariant space V (see section l2~2j) . The 
genus of the curve r is g, where g is as in (fTTjl . 

Proposition 3.3. Let k = max(Z , h, h, h)> ^3 = m i n (^o ; hi h) an d I — lo + h + 
h + h- If Conjecture^ is true, then the arithmetic genus of the curve T is given by 

ko, I is even and ko + k 3 > 1/2; 

(/ — 2k 3 )/2, I is even and ko + k 3 < 1/2; 

k , I is odd and k > (I + l)/2; 

(/ + l)/2, I is odd and k < (l + 1) /2. 

Proof. It follows directly from Proposition 12.31 □ 

Remark. The values in (|T%|) have already appeared in jlJ^JCSl, although the present 
author could not understand the proofs of computing the genus of V given in |H IT4"] 
because he could not confirm how to avoid the possibility that either of the equations 
P{E) = and Q{E) = could have multiple roots. 

By combining Propositions 12.61 and 13.31 we obtain the following corollary. 

Corollary 3.4. If two of li (i = 0,1,2,3,) are zero, then the arithmetic genus of the 
curve T is max(Z , h, h, h)- 

Proposition 3.5. If Iq = l\ — 0, I2, h £ Z>o an d r G Z + y/— lR>o, then u(x) = 
Y^i=2 h{h + l)p( a; + ^i) is a finite-gap potential. 

Proof. First from the condition r G Z + \J — 1K. >0 we have u{x) G R for x G R. 
Also, since the poles of the function p(x + r/2) (resp. p(x + (r + l)/2)) are located 
on the set Z + (Z + l/2)r (resp. (Z + 1/2) + (Z + l/2)r ), the function u(x) is 
holomorphic on R. Finally it is shown in fTJ |2] that, if u(x) is real-holomorphic on 
R and u(x + 1) = u(x), then u(x) is a finite-gap potential if and only if u(x) is an 
algebro-geometric finite-gap potential. Thus we obtain the proposition. □ 
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3.2. The determinantal formula for the commuting operator. 

In this subsection a determinantal formula for the commuting operator A is intro- 
duced. 

Conjecture 2. Let U (i = 0,1,2,3,) be nonnegative integers and V be the invariant 
subspace defined in section WJ\ 

Let fi(x), . . . , f2g+i(x) be a basis for the space V . Then 

( fl(x) ... f2g+l(x) 1 \ 

(£) AO*) ••• (£) fmiW (£) 



(19) 



A = Ar 



V(±)*"/i(*) 



, dx J 



J 



for some nonzero constant A Q . Here the determinant for the matrix 



• • • a l,n 



Qj n \ . . . Q"n,n 



is defined by J2 ae&n sgn(a) 

Proposition 3.6. If two of Zj (i = 0, 1,2,3) are zero, Conjecture^ is true. 

Proof. It is sufficient to show the statements for the case l 2 = h = and l > l\. 
From Propositions 12 . 61 and 13 . 21 we have A 2 = —P(H). In Theorem 3.2] it is shown 
that the equation P(E) = has no multiple roots for generic r. For this case we 
have dim\^ = 2g + 1, where g is defined in Proposition 12.11 

We fix the period r to satisfy the condition that the equation P(E) = has no 
multiple roots. Let E{ (i — 1, . . . , 2<?+l) be the zeros of P(E) = and let Vi be the 2- 
dimensional space spanned by solutions to the differential equation (H — Ei)f(x) = 0. 

From the commutativity [A, H] = the operator A preserves the space V%. For 
f(x) G Vi we have A 2 f(x) = -P(H)f(x) = -P{E i )f{x) = 0. Therefore all eigenval- 
ues of A on are 0. 

Let £i(x) G V H Vi. From the definition of V the function £i(x) is doubly periodic 
up to signs and satisfies the equation (H — E^^x) = 0. Another solution to (H — 
Ei)f(x) = is given by £i(x) J Q X dy / (,i(y) 2 ■ This function could never be a doubly 
periodic function up to signs because it has Iq + 1 zeros at x = and has fewer poles 
on the fundamental domain. Hence the dimension of the space of doubly periodic 
functions up to signs satisfying the differential equation (H—Ei)f(x) = is one. Since 
the function E) is doubly periodic, the coefficients of (d/dx) 3 (j — 0, . . . , 2g + 1) 
on A are doubly periodic (see (|T2"j)): hence the function A£i(x) is proportional to 
&(x). Thus we have A£i(x) = 0, because the eigenvalue must be zero and therefore 
we also have Af(x) = for all f(x) G V. 

Since degA = dimV, all solutions to Af(x) = must satisfy f(x) G V and we 
have 

tM*) ■■■ 



(20) 



A = h(x) 



( 



1 
d_ 

dx : 

(AY 

\dx) 



\ 



x 



(If 3+1 h^) (if* 1 ) 
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for some function h(x). Since the coefficient of (d/dx) 2g+1 on the right-hand side of 
(J201) is (-1) 9 (see l)TZjl ). we have 



(21) 



/ hi*) 



d h(x) 



hg+l{ 



dx 



x 



_x) 

£f2g+l(x) 



\ 



h(x) 



By differentiating equality ([2*T]). we obtain that the coefficient of (d/dx) 29 on the right- 
hand side of (|2U|) is equal to (—l) 9 h'(x)/h(x). From the commutativity [A, H] = 
the coefficient of {d/dx) 29 in the operator A must be zero, whence h'(x) = 0. 

Therefore relation (|19J) follows for the case 1% = Z3 = and Iq > h, if the equation 
P(E) = has no multiple roots. From analytic continuation in r we obtain relation 
(Unj) for all r e R + v / ^lR>o when l 2 = h = 0, Z > ii- □ 



3.3. Monodromy. 

In this subsection we calculate the monodromy, specifically the constant B(E) 
for A(x + 1,E) = B(E)A(x, E), where the function A(x,E) is the solution to the 
equation (H — E)A(x,E) = of the form The monodromy plays an important 
role in investigating eigenvalues and eigenstates for systems with physical boundary 
conditions. 

From formula (EI) we have 



(22) 



rl+e 

A(x + 1, E) = ±A(x, E) exp / 

JO+e 



-Q(E)dx 



with e a constant so determined as to avoid passing through the poles in the in- 
tegration. The sign ± is determined by the analytic continuation of the function 

v /!(Z7a 

Although this formula is explicit, the calculation of the integral would be difficult. 
In [H)] it is shown that the function A(x, E) admits the following expression of the 
Bethe ansatz form: 



(23) 



A(x,E) 



exp(7ri/— Tex) Y\ l j=i @( x + tj) 



0{x) l »6(x + l/2) l ^9(x + (1 + r) /2) l *0{x + t/2) 1 * 



where I = lo + h + h + h and 9(x) is the elliptic theta function. The values tj 
(j = 1, . . . , I) and c satisfy the Bethe ansatz equation given in ^f)J Theorem 3.13]. 
From this expression we have 



(24) 



A(x + 1, E) = A(x, E) exp(7r\/ z Tc). 



To obtain the value c we must solve the Bethe ansatz equation, which is in general 
transcendental. 

In this subsection we introduce a different formula. The monodromy is expressed as 
an hyperelliptic integral of second kind. The following formula for a smooth finite-gap 
potential is found in P ITUj. 



Theorem 3.7. Let k e Z> 0; P G 2Zu 1 ®2Zuj 3 (= Z©Zr) 
of H such that A(x + P, E { Q P) ) = A(x, E { P) ) (resp. A(x + P, E ( P) ) 



( P) 

and Eq be the eigenvalue 



-A(*,4 P) )) 
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Then 

(25) A(x + P, E.) = A(x, E.) exp ( -- J dE 

, 1 f E * L P + £ Z(x,E)dx 

resp. A(x + P, E m ) = -A(x, E*) exp — / Jo+£ —dE 



2 J E p y^Q(E) 

with e a constant so determined as to avoid passing through the poles in the integra- 
tion. In particular, if E is the eigenvalue of H such that A(x + 1,Eq) = A(x,E ) 
(resp. A(x + l,E ) = —A(x,E )), then 

(26) K{x + l 1 E m )=A(x 1 E m )exp(~ f" ^L dE 

\ l JEa \/-Q{E) 



E y/-Q(E) 

resp. A{x + 1, E.) = -A{x, E.) exp ^~ ^ ^^_ dE 
where the polynomial Qi(E) of degree g is given by 

rl+e 

(27) Qi(E)= / E(x,E)dx. 

JO+e 

Proof. Set 

(28) R = R(x,E) = ^^L 

y/—Q{E) 

We use ' instead of d/dx. Set v = J2i=o h(k + ^)p{ x + ^i)- From (j5J) 

(29) R 2 (v — E) — ^RR" + ~ (R'f = 1. 

By differentiating (f2H|) with respect to E we have 

/ x dR" R'dR' nn R" OR Al „.dR 

Then 

. . / dR' R'dRY 
(31) {-M + R8e) 



R dE R dE v ' dE) \R R 2 J dE R dE 

= 2R+ R L ^ ( 2RR " ~ AR2{V ~ E) - {R,)2) 

= 2R-± 2 dR =2R + 4JL(± 
R 2 dE dE \R 

Since the function — + ^f§ is doubly periodic in x, it follows that 



f P+£ ( dR ' R'dR\' n 
(32) / ( -— + — — ) dx = 0. 

JO+e 



dE R dE 
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Hence 

i r p+£ r p+£ d /i\ 

(33) -2L Rdx rL m UK 



By integrating with respect to £/ we see that 
(34) exp 



( rP+e y/~Q(E), \ 



From formula (jBJ) we have 
(35) A(x + P,E) = ±A(x, E) exp ^Z?^!^ ■ 



The sign ± is determined by the analytic continuation of the function a/2(x, E). By 
the assumption A(x + P,E { P) ) = A(x,E (P) ) (resp. A{x + P,E (P) ) = -A(x,E {P) )) 



we have ±exp J 0+£ = 1 ( res P- ±ex P Jo+e S ( x ,e^) = ^ W ere the 

sign ± is chosen to be the same as the one in ()35|) . Combining with (jMj) and (j33j) we 
obtain □ 

Remark, (i) A similar equality to (}3*Tj) can be found in [3]. 

(ii) Although the integral J^* -Q^SL= dE may be dependent on the choice of the initial 

point Eq such that A{x + 1, E ) = A(x, E ) (resp. A(x + 1, E ) = —A(x, E )) and the 

path from E to E*. the function exp ( — \ fif* 'j 1 ^ dE ) +exp ( ~ f,f* -^2=dE ) 

V. y/-Q(E) J V y/-Q{E) J 

is independent on £7 with A(x + 1, £"o) = A(x, £"o) (resp. A{x + 1, £? ) = — A(x, E )) 
and the path. 

From Proposition 123 if Q{Eq) = 0, then the differential equation (H — Eo)f(x) = 
(see (J2J)) has a solution of doubly periodic function up to signs. It is seen that the 
solution is written as A (a;, Eq) = ^/S(x, Eq) up to constant multiplication. Hence 
A(x + 1, E ) = A(x, E ) or A(x + 1, £ ) = -A(x, Eq). Combined with Theorem E3 
we obtain the following corollary: 

Corollary 3.8. Let l{ G Z>o and let Eq be a solution to Q(E) = 0. Then A(x + 
1,E ) = A(x,E ) or A(x + 1,E ) = -A(x,E ). IfA(x + l,E ) = A(x,E ), then 

A(x + 1, E m ) = A{x, E.) exp (-\ jj; J^LdEj and, if A(x + 1, Eq) = -A(x, Eq), 
then A(x + 1, E.) = -A(x, E.) exp (-1 jj ^g^) ■ 

Theoretically the integral Qi(E) = Z(x, E)dx can be calculated explicitly. 
This is illustrated by some examples in the next subsection. 

3.4. Examples. 

Two examples are provided to illustrate plainly what has been done to this point. 
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3.4.1. The case l = 0, l\ = 0, l 2 = 0, l 3 
The Hamiltonian H is given by 



KOUICHI TAKEMURA 
1. 



(36) 



d 2 ( T 



For this case the space of finite-dimensional elliptic functions V defined in section E21 
is spanned by co-p functions p\(x + |), p 2 (x + |) and p 3 (x + |). Hence dim = 3. 
The characteristic polynomial of the Hamiltonian H on the space V is 



(37) 



P{E) = (E + ei ){E + e 2 )(E + e 3 ), 



where p{\) = ei, p{^y~) = e 2 , and p(|) = e3. For this case we have Q(E) = P(E) 
from Proposition 13.61 

From Proposition ^. 5l 2ofx+ 1) is a finite-gap potential if r G Z- 



. Indeed, 



if r G 2Z + 
(38) 



-HR>o, then we have 



\ a b (H) = (-oo, -ex) U (-e 2 , -e 3 ) 



for -ei < -e 2 < — e 3 (see Definition [TJ). If r G 2Z + 1 + 
— e 2 and 



then — ei < — e 3 < 



(39) R\a b (H) = (-oo,-e 1 )U(-e 3 ,-e 2 ). 

The functions -E 1 ) and A(x, i?) are given by 



(40) ~(x,E) = p(x+^+E, 
A(x,E) = Up (x + ^j 



+ E exp 



x ^-(E + ei )(E + e 2 )(E + e:i ) 



dx. 



p(x + |) + £ 

and -E) satisfies 



Here A(±x, E) satisfy the equations (H — E)A(±x, E) 
equation (j3J). 

The function 2p(x + |) is also an algebro-geometric potential if r G R + 
The commuting operator A (|T2*j) is given by 

3 



■1R 



>o- 



(41) 



A 




3p I x + 



d 
dx 



As in the case of Proposition 13.61 we obtain the following determinantal formula: 

(42) (ei-e 2 )(e 2 -e 3 )(ei-e 3 )A 

/ Pl( X + l) P2(X + ^) p 3 (x+l) 1 \ 

p[(x + ^) p' 2 {x + l) p' 3 (x + $ £ 

+ Pg(x + i) p?(x + 5) (£) 2 

For the operators A and if there is an algebraic relation 

(43) A 2 = -(H + ei )(H + e 2 )(H + e 3 ) 
and the genus of this curve is one. 
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Now we can proceed to calculate the coefficient of the monodromy We have 



(44) / E(x,E)dx = [ (E + p(x + ^-))dx 

= E-((l + ^ + e)+((^ + e) =E-2 Vl , 

where r\\ = C(|)- (Here we used definitions and formulas from the appendix.) Since 
E = — ei is a solution to Q(E) = such that A(x + 1, E ) = A(x, E ), it follows that 

(45) L(x + l,E m ) = A(x,E,)e*p(~ dE) , 

V 2 J-ei y / -(E + e 1 )(E + e 2 )(E + e 3 ) J 

which is an elliptic integral of second kind. 



3.4.2. The case l = 2, l x = 1, l 2 = 0, l 3 = 0. 
The Hamiltonian H is given by 



(46) H = + 6p(x) + 2p fx + ~ 



For this case the space of finite-dimensional elliptic functions V defined in section 
12. 21 is spanned by the functions p(x) — e x , p 2 (x)/pi(x), p 3 (x)/pi(x), pi(x)p 2 (x), and 
pi(x)ps(x), so dimV^ = 5. The characteristic polynomial of the Hamiltonian H on 
the space V is 

(47) P(E) = —(E + 4 ei ){16E 4 - 224 ei E 3 + (56g 2 - \Ue\)E 2 

16 

- (200# 2 ei - 5920e?)£ + (81^ - 3640# 2 e? + 19216e?)}, 



where g 2 = — 4(eie 2 + e 2 e 3 + e^ei). For this case we have Q(E) = P(E) from 
Proposition 13.61 

The functions S(x, E) and A(x, E) are given by 



(48) 

S(a;, E) =9p(x) 2 + 3(E - 2e x )p{x) + (E + 4ei)p (x + + E 2 - he x E - 21e\, 

A(x, E) =y/E(x, E)exp f ^— ®^ dx. 

J z.[x,E) 



Here A(±x,E) satisfy the equations (H — E)A(±x,E) = and S(x,E) satisfies 
equation (jSJ). 
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The function 6p(x) + 2p(x + |) is an algebro-geometric finite-gap potential. The 
commuting operator A (|12|) is given by 



d \ 5 _ / , ( 1\ \ ( d^ ' 




< 49 » A =UJ - 5 r (l)+p l l+ 2j" c vv* 
_»( vw+ ^(. + »))(-)• 

- 15 Up(x) 2 + 2p(x + ^ - 2e lP (x + - 3e? - ^ 2 

- 15 fx + ij - p' fx + i 

For the operators A and if there is an algebraic relation A 2 + Q(H) = and the 
genus of this curve is 2. 

As in the previous example we calculate the coefficient of the monodromy (J26j) . 
From the formulas p(x) 2 = ^p"(x) + C'( x ) = ~pi. x )i C( x + 1) ~ C( x ) = ^Vi an< l 
the periodicity of p'(x), we have 

rl+e 3 

(50) / E(x, E)dx = E 2 — (5ei + 8 Vl )E + -g 2 - 27ej + Ae^. 

Since E = — 4ei is a solution to Q(E) = such that A (re + l,E ) = A(x,E ), it 

follows that 

(51) 

\ 2 "'-4ei y/—Q{E) 

which is an hyperelliptic integral related to the genus two curve. 

4. Application to the spectral problem 

4.1. The strategy. In this section eigenvalues of the Hamiltonian H for the systems 
with physical boundary conditions are considered. More precisely we consider eigen- 
states and eigenvalues of the Hamiltonian in an appropriate Hilbert space. In the 
general setting the Hilbert space is taken as the space of functions that are square- 
integrable for some inner product. For the infinite-dimensional case the Hamiltonian 
generally does not act on all elements of the whole Hilbert space but only on a dense 
subspace. Hence we should specify the Hilbert space and its dense subspace for each 
model. 

In section 14.21 the Hilbert space and its dense subspace for the BC\ Inozemtsev 
model are introduced. In previous sections we obtained the eigenfunctions A(x, E) 
of the Hamiltonian if in © without boundary conditions and the constant for the 
quasi-periodicity in (|26|). It is mentioned in this section how the conditions for square- 
integrability correspond to the ones for quasi-periodicity. 

For the trigonometric case, i.e. r — > a/— loo, the model becomes much simpler. In 
this case all eigenvalues and eigenfunctions are described explicitly. The eigenvalues 
of the nontrigonometric models are obtained by connecting from the trigonometric 
model. Here the Hamiltonian H depends upon the variable p = exp^V— It) and 
the trigonometric case corresponds to the case p = 0. It is shown in that the 
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eigenvalues are holomorphic in p if — 1 < p < 1 by an application of the theory of 
holomorphic perturbation due to Kato and Rellich jU] . 

To further understand the eigenvalues we examine the singular points of the eigen- 
values in the variable p on the complex domain \p\ < 1. In section l4~3*l a sufficient 
condition for the holomorphy of an eigenvalue at a given point p = po is provided. 
Here relation ()26|) plays an important role. 

Firstly the case l G Z>i and l\ G Z>i is discussed in sections 14.21 and 14.31 In 
section 14.41 we treat the case Iq G Z>i and l\ = 0, although the case Iq — and 
/i G Z>i can be treated similarly. In section |4~B1 the case Iq = 1% = is discussed. 

4.2. The Hilbert space and the trigonometric limit. 

In this subsection assume that l , li G Z>i and / 2 , Z3 G Z> and set p = exp(7r-\/— Tr). 
Then the potential of the Hamiltonian H is dependent on p. We introduce the Hilbert 
space and its dense subspace. Note that our notation differs from that in ^7]. Define 
the following function, which is indeed the ground state for the trigonometric (p = 0) 
model: 

(52) $(ac) = (sin7rx) fa+1 (cos7rx) il+1 . 

By using the function $(x), we define the Hilbert space H and its subspace D for 
every p (\p\ < 1) as follows: 

£ \f( x )\ 2dx < +°°> 
f( x ) = (-l) l °+ l if(x + l) a.e. x, 
f \ x ) = (— l) i o+i/(_ x ) a .e. x 



(53) H = { f : R -»• C measurable 
and 



(54) D= /: 



/(») 

/(2) = (-1)^/^ + 1), /(z) = (_l)'o + l /( _ x) 



is holomorphic on 



Observe that the spaces H and D are independent of p. It is shown that the space 
D is a dense subspace of H. 

Remark. Let 

J 1 |/(x)| 2 |$(x)| 2 ^<+oo, 



H ' CUU ' aMll ' a /(x) = /(x + l), f{x) = /(-x) a.e. x 

be the gauge-transformed Hilbert space. (This space is denoted as H in |17j.) There is 
a simple relationship between the spaces H 5 and H, namely h(x) G H 9 <^ /i(a;) ( l ) (x) G 
H. 

It is shown in that for each eigenstate ip m (x) G H of the Hamiltonian H we 
have ip m {x) G D if p G R and |p| is sufficiently small. Thus it would be reasonable to 
consider the spectral problem for the space D. 

Proposition 4.1. The Hamiltonian H preserves the space D for every p G C such 
that \p\ < 1. 

Proof. Let p be a complex number such that |p| < 1 and let TC(p) = $(x) _1 oifo$(x) 
be the gauge-transformed Hamiltonian. Then we have 

. , . , d 2 N cos7rx N sin7rx\ . , 
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where V p (x) is an even periodic function in x that is holomorphic onieR (for detais 
see [T7J §2]). It is enough to show that the operator T~C(p) preserves the space D g , 
where 



(56) r> g = \f:R^C 



f(x) is holomorphic on R, 
f(x) = f(x + l),f(x) = f(-x) 

Let f(x) G D. Then the function f(x) = TC(j))f(x) satisfies f(x + 1) = f(x) = 
f(—x) and f(x — 1/2) = /(1/2 — x). The function f(x) is meromorphic on R and 
may have poles at x G |Z of degree at most 1. If /(x) has a pole at x = (resp. 
x = 1/2), it contradicts the symmetry f(x) = f(—x) (resp. f(x— 1/2) = /(1/2 — x)). 
Hence the function f(x) has no poles at x = and x = 1/2. Furthermore from the 
periodicity the function f(x) has no poles on R thus completing the proof. □ 

For the case p = the problem becomes much easier. We briefly describe the 
eigenvalues (for details, see [T7|). 

If p — > (r — > -\/ — I 00 )) then we have p(x) — > — 7T 2 /3 + 7r 2 / sin 2 7rx etc. and 
H-^H T -n 2 (EL W + !)) A where 

,j2 2 2 

(57) # T= +z (z + i)_Z_ + /l (/ 1 + r "" 



<3X Z Sin 7TX COS z TTX 

The (m + l)st lowest eigenvalue of the trigonometric Hamiltonian H? on the space 
H is 7r 2 (2m + Zq + ^i + 2) 2 . The corresponding eigenfunction i/j m (x) is written as the 
product of the ground state $(x) from (J52*)l and the Jacobi polynomial of degree m. 
In particular we have ip m (x) G D. 

If — 1 < p < 1, then the existence, completeness, and real-holomorphy of the 
eigenvalues in p are shown in Theorem 3.6] by application of Kato-Rellich theory. 

Proposition 4.2. (c.f. [T71 Theorem 3.6]) Let l ,li G Z>i and h-,h G Z> . ^4// 
eigenvalues of H for — 1 < p < 1 on £/ie Hilbert space H can 6e represented as 
E m (p)(= 7r 2 (2m+/ + /i + 2) 2 — ^ Z}i=o^(^+l)) ( m e ^>o); w/wc/i real-holomorphic 
in p G (—1, 1), and i/ie uaZue E m (0) coincides with the (m + l)st lowest eigenvalue of 
the trigonometric model. 

Although the eigenvalue E m (p) is real-holomorphic, it may have singularities on 
the complex domain. This is investigated in the next subsection. 

4.3. The condition for square-integrability. 

In this subsection we write down the condition to determine which values of E 
admit square-integrable eigenf unctions and discuss the singularities of the eigenvalue 
E m {p) in the variable p = exp(7T\/— Tr). Throughout this section it is assumed that 
l , h G Z>! and l 2 , l 3 G Z> . 

Let A(x, E) be the function defined in (JBJ). Then HA(x,E) = EA(x,E) and 

A(x+ 1,E) = B(E)A(x,E), where B(E) = exp (-± -Qj^dEJ, E is the 

value satisfying A(x + 1,E ) = A(x,E ) and Q\(E) = f Q ^ H(x, E)dx (see Theorem 
E2|). 

Set A s y m (x,E) = A(x,E) - (— 1)'°A(— x, E). Then the relation if A sym (x, E) = 
EA sym (x, E) is obvious. The following proposition gives a sufficient condition for 
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the square-integrability because, if A aym (x, E) G D, then the function A aym (x, E) is 
square-integrable on (0, 1). 

Proposition 4.3. (%) If B(E) = (-l) lo+h , then A^ m (x,£) G D. 
(ii) If Q(E) ^ 0, then the function A sym (x,E) is nonzero. 

Proof, (i) We must check that the function A sym (x, E) satisfies the defining relations 
(PI) for D. 

The relation A sym (x,E) = (— iy° +1 A sym (— x, E) is trivial from the definition of 
A sym (x,£). The relation A sym (x,£) = (-l)'o+'i A sym (x + 1,E) is obtained from the 
condition B(E) = (— l) io+il . Since the function A sym (x, E) satisfies the differential 
equation (J2J), it may have poles only at x G \L on R. 

We have the expansion A sym (x, E) = x a (c + CiX + ■ ■ • ), where c ^ and (a = — 1 
or l + 1), because the exponents of the differential equation (j2J) at x = are — 1 and 
Z + 1. From the property A sym (x,£) = (— l) /[,+1 A sym (— x, E) and Z e ^>o, we have 
a = /o + 1. Thus the function A sym (x, £7) is holomorphic at x = 0. 

We set A(x) = A sym (x - 1/2, E). Then we have A(x) = A sym (x - 1/2, E) = 
(-l)'°+ Zl A sym (x + 1/2, £) = (-l) /l+1 A sym (-x - 1/2, E) = (-l)' 1+1 A(-x). Similarly 
to the case A sym (x, E)\ X=Q , the function X(x) is holomorphic at x = 0. Hence the 
function A sym (x, E) is holomorphic at x = 1/2. 

From the periodicity of A sym (x, E) we obtain that the function A sym (x, E) is holo- 
morphic on R. Therefore we have checked all the defining relations and have proved 
«. 

(ii) It is sufficient to show that the functions A(x, E) and A(— x, E) are linearly inde- 
pendent. This is obvious from the expressions A(— x, E) = a/5(x, E) exp J 



Proposition 4.4. Let l , l x G Z> x and Z 2 , Z 3 G Z> . //5(E) ^ and Q(£) ^ 

0, then any nonzero solution to (0) is not square-integrable. In particular any value 
E such that B(E) ^ (— l)' 0+il and Q{E) ^ is not an eigenvalue of the Hamiltonian 
H in the space D. 

Proof. From the condition Q(E) ^ any solution to (J2J) can be written as a linear 
combination of A(x,E) and A(— x, E). From the inequality B(E) ^ (— l) Zo+?1 it is 
impossible to remove the poles at x — 0, 1/2, and 1 simultaneously. □ 

For the case Q(-E) = a different argument is required. 

We observe the condition under which eigenvalues in D are analytically continued. 
The eigenvalues in D are determined by the condition B(E) = (— l)' 0+il . The function 
B(E) depends on two variables, E and p, and furthermore is holomorphic in E and 
p provided (E,p) G C x B x and Q(E) ^ 0, where B x = {p G C | \p\ < 1}. 

Let \p#\ < 1 and suppose that B(E*) = (— l)' 0+il at p = p*. By the implicit function 
theorem, if ^\ p=Pt E=Et ^ and Q(E*) ^ 0, then there exist e G M>o and a unique 
holomorphic function E(p) on {p G C | |p — p*| < e} such that B(E{p)) = (— l)'°+'i 
and E(pt) = E*. Thus we have the following proposition. 

Proposition 4.5. Lei |p*| < 1 and fre an eigenvalue of H on D snca t/iat 
§§ Ip=p,,.e=e, 7^ and Q(E*) ^ 0. Tnen i/iis eigenvalue on D zs analytically con- 
tinued near p = p* and it is holomorphic at p = p*. 




□ 
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From formula ()26|) we have more explicit results. 
Theorem 4.6. Let lo,h G Z>i and h,h G Z>o- 

(58) / * -Q^2=dE G 2(/ + /Ottv/^T + 47T V /Z TZ, 

the value E satisfies A(x + l,E ) = A(x,E ), Qi(E) is a polynomial in E defined 
by Qx{E) = J^E(x,E)dx and Q(EJ ^ 0, then we have A s y m (x,E,) G D and 

A s y m (x,E*) o. £ 

(ii) If \p*\ < 1 and E* is such that Qi(E^) ^ 0, Q(E*) ^ and there exists a nonzero 
function f(x) G D satisfying (H — E*)f(x) = at p = p*, then the eigenvalue E* is 
analytically continued near p = p* and is holomorphic at p — p*. 

Proof. From (|2T)j) the condition B(E*) = (— l)' ()+/l is equivalent to the condition (j58j) . 
In addition, as is shown in Propositions 14.31 and 14.41 if Q{E*) ^ 0, then the condition 
that there exists a nonzero function f(x) G D such that (H — E*)f(x) = is also 
equivalent to B(E+) = (-l) lo+h . 

Finally it is easy to see that the conditions ^\ p=Pt ^ E=Et ^ and Q(E*) ^ are 
equivalent to the conditions Qi(E*) ^ and Q(E*) ^ 0. By combining these facts 
and Proposition 14.51 we obtain the theorem. □ 

Remark. If the value E satisfies A(x + 1, E ) = —A(x, E Q ), then the condition (JSHJ) 
is replaced by 

(59) / -Q^E)=dE G 2(Z + h + + 4ttv /Z TZ. 
Je* ^ 



>e ^q[e) 

Let E m (p) be the eigenvalue of the Hamiltonian H on the Hilbert space H defined 
in Proposition 14.21 Then from [T7j the eigenfunction v m (x,p) of the eigenvalue E m (p) 
on the full Hilbert space H satisfies v m (x,p) G D, if p G M and \p\ is sufficiently small. 

Hence jj m(p) -^=L=dE G 2(Z + h)^V^l + Any^lZ, if E satisfies A(x + 1, E ) = 

A(x,Eq), p G K and |p| is sufficiently small. From Proposition 14.21 the eigenvalue 
E m (p) is holomorphic if —1 < p < 1. 

From Theorem 14.61 if Q(E m (p)) ^ and £7 satisfies A(x + 1, E ) = A(x, E ), then 
the eigenvalue E m (p) is so continued analytically in p(e C) as to maintain the relation 
rEmip) Qi(E) dE 2 n + / 1 ) 7rv ^T + 47T V / ^TZ. Let p m G C such that IpJ < 1. If 

JE o y/-Q{E) y ' * ^ 1 

Q(E m {p*)) 7^ and Qi(-E m (p*)) 7^ 0, then the eigenvalue E m (p) does not have a 
singularity at p = p*; in other words the eigenvalue E m (p) can have a singularity at 
p = p* only when Q(E m (p*)) = or Qi(Ejj) m )) =0. 

As an aside note that in an algorithm for calculating the eigenvalue E m (p) 
as a power series in p was introduced and the convergence of the power series for 
sufficiently small \p\ was shown. Numerically it can be demonstrated that the radius 
of convergence of the power series is strictly less than 1. Assume that the radius of 
convergence of the power series is p c G (0,1). Then from the preceding discussion 
there exists a solution to Q(E m (p*)) = or Q\(E m {p*)) = satisfying |p*| = p c . In 
this way the condition (Q(E) = or Qi(E) = 0) is closely connected to the radius 
of convergence of the eigenvalues. 
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For the case Iq — 2, li — 1, l 2 — and Z 3 = the polynomial Q(E) is given 
in (|47|) and the polynomial Qi(E) is given as the right-hand side of (J5(J|) . Then 
we have deg E Q(E) = 5 and deg E Qi(E) = 2. To find the singular points of the 
eigenvalues for the space D it would be effective to compute the values E* such that 
(Q(E*) = or Qi(E*) = 0) as varies and investigate when these values E* also 
satisfy f_J ei J l( ^dE E (4Z + 2)tt v / ^T. It is anticipated that more such results will 

be reportable in the future. 

4.4. The case l E Z>i and h = 0. 

In this subsection we consider the case Iq E Z>i, l\ = 0, h, h £ ^>o- For this case 
consider the larger Hilbert space H* and its dense subspace D* defined by 

f*\f(x)\ 2 dx < +oo, 
fix) = fix + 2) a.e. x, 
f(x) = (_iyo+ij(_ x ) a e x 



(60) 



C measurable 



and 



(61) 



D* 



7 ■ ^L+i is holomorphic on 

f(x) = f(x + 2), /(x) = (- 



-V(- 



-x) 



Then D* C H*, H C H*, and D C D*. It can be shown that the Hamiltonian H 
preserves the space D* (c.f. Proposition 14. lj) . 

For the case l E Z> 1; Z x = 0, l 2 , ^3 £ %>>o, & U eigenvalues of if for —1 < p < 1 on the 
Hilbert space H* can be represented as E^Jj)) (m E Z> ), which is real-holomorphic 
in p E (-1, 1) and for which ££(0)(= vr 2 (m + Z + l) 2 - ^ Ei=o W + 1)) coincides 
with the (m + l)st lowest eigenvalue of the trigonometric model (see ^3 Theorem 
3.6]). 

Let A(x, E) be the function defined in (JHJ). 



A(x + 1,£) = B(E)A(x,E), where 5(E) 



Then HAix, E) ■■ 
e^i-ijfJMH 



E 

E ° y/-Q(E) 



EA(x,E) and 
dE | , Eq satis- 



fies A(x + 1,Eq) = A(x,E ), and Qi(E) = J^E(x,E)dx. Set A s y m (x,£) = 
A(x,E) - (-l) ,0 A(-ar,S). If B(E) = 1 or -1, then we have A s y m (x,£) G D* and, 
if Q(E) 7^ 0, then the function A sym (x, E) is nonzero. Conversely, if B(E) ^ ±1 and 
Q(-E') 7^ 0, then it can be shown that any nonzero solution to (J2J) is not square- 
integrable. Note that the condition B[E) = ±1 is equivalent to the condition 
rE q^e) d £ e 2ny /^iz, where E satisfies Q(E ) = 0. 

Similar statements to Proposition 14.51 and Theorem 14.61 can also be derived. 



4.5. The case l = l x = 0. 

In this subsection we consider the case Z 
Set 



(62) 

and 
(63) 

Then D** 



h* 



measurable 



h — 0, I2, h £ Z> . 



£ \f( x )\ 2dx < + °°> 
/(x)=/(z + 2) a.e. 



D 



C H* 



/:R- 
Hcffcff 



fix) is holomorphic on 

f{x) = f{x + 2) 

and DcD*C D**. 
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It is shown that the Hamiltonian H preserves the space D** (c.f. Proposition 14. lj) . 

For the case Iq — l\ — 0, l 2 ,h £ ^>o, all eigenvalues of H for — 1 < p < 1 on 
the Hilbert space H** can be represented as E%(p) (m G Z> ) and E^(p) (m G 
Z>i), which are real-holomorphic in p G (—1,1) and for which -E^(O) = -n 2 m 2 — 

jY^=2 l *{k + 1) (m G Z> ) and £°(0) = vr 2 m 2 - f £^2 W + 1) (™ e Z>i) (see 
[T7| Theorem 3.6]). 

Let A(x, E) be the function defined in (jSJ. Then HA(x, E) = EA(x, E) and A(sc + 
1,£) = B(E)A(x,E), where 5(E) = exp -^SLdEJ, E satisfies A(ac + 

1, £/o) = A(x, Eq) and Qi(E) = H(x, E)dx. It is easy to show that, if 5(E) = ±1, 
then A(x, E) G D**. The condition B(E*) = ±1 is equivalent to the condition 

For the case Iq — l\ — 0, l 2 , h £ ^>o we have the following proposition. Some parts 
are proved similarly to Proposition 14.51 and other parts are proved directly. 



Proposition 4.7. (%) If Q(E*) ^ and B(E+) = ±1, then the functions A(x,E*) 
and A(—x, E*) are linearly independent and thus E* is a double eigenvalue, 
(ii) If Q(E*) = 0, then the function A(x,E*) is doubly periodic up to signs and the 
relations (A(x,E*) = A(—x,E*) or —A(—x,E*)) and B(E*) = ±1 are obtained. In 
this case the eigenfunction A(x,E*) corresponds to the nonrepeated eigenvalue E*. 
(Hi) The nonrepeated eigenvalues are obtained by the diagonalization of the finite- 
dimensional space V defined in section WM 

(iv) The repeated eigenvalues E* are analytically so continued in p as to preserve the 
property ff* J l{E) dE G 2n^lZ, where E satisfies Q{E ) = 0. // Qi(E*) + 

and Q(E*) ^ at p = p*, then the eigenvalue is analytic in p at p = p*. 



By application of Proposition 14. 71 the following statements related to the finite-gap 
property can be shown: 

• If l 2 + h is even and — 1 < p < 1, then E 2m {p) = E 2m {p) f° r m > §(^2 + h) an d 

Efm-iiP) = ^2m-i(p) f° r m > W 2 ~ ^ 3 I' ^ ne res ^ °^ ^ ne eigenvalues correspond to 
the ones of the Hamiltonian H on the space V defined in section 12.21 

• If I2 + h is odd and — 1 < p < 1, then E^Jj)) = E2 m (p) for m > h(\h — h\ — 1) and 
Etm-xip) = E2 m -i(p) for in > ~(l 2 + ^3 + 1)- The rest of the eigenvalues correspond 
to the ones of the Hamiltonian H on the space V defined in section 12.21 

For the case Iq = li = l 2 = and I3 = 1 the results are illustrated more ex- 
plicitly. In this case Q(E) = (E + e x ){E + e 2 )(E + e 3 ) and Qi(E) = E - 2^. 
From Proposition 14.71 the nonrepeated eigenvalues are — e\, — e 2 and —e^. Moreover 
Eq(p) = ~ e i an d {E>i(p), Ei{p)} = {-62,-63} whereas, for the other eigenvalues, 
Emip) — E^{p) (m G Z> 2 ). Finally the double eigenvalue E°(p) (m G Z> 2 ) is ana- 
lytically so continued as to preserve f Em ^ ^ =dE g ItxJ —YL and, 

J _ J 1 J ' e ^ ^/-{E+e 1 )(E+e 2 ){E+e i ) V 

if E^(p*) 2771, -ei, — e 2 or -e 3 , then E°(p) is analytic at p — p*. 
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5. The Heun equation and the Inozemtsev model 

In this subsection the relationship between the Heun equation and the Inozemtsev 
model is reviewed and the results in this paper are translated to the ones for the 
Heun equation. 

Let f(x) be an eigenfunction of the Hamiltonian H with the eigenvalue E (j2J), i.e., 



3 



dx 2 



i=0 



e 2 - e 3 



$(u>) 



(64) 
Set 

(65) a 

ei - e 3 ' 
Under the transformation 

(66) " - pM 

equation (JMj) is equivalent to the equation 
(67) 

where /(^^KI><^^ 



+ > l i (l i + l)p(x + cj i )-E\f(x) = 0. 



w 



w 2 (w 



ei - e 3 



I) * (aw 



C3 




d 



+ 



af3w — q 



f(w) = 0, 



- p(x)-e 3 I * V p( x )-e 3 

h + h-h + 3), 7 = + i <) 



-, a 



h + 



dw ' w(w — l)(w — t) 

~(h + h + h + h+ ^), (3 = ±(Z + 

+ c ) , and 



l 2 + lv = -i 



2> "ii2' • 2 1 ~* 4o \ ei-e 3 

Co = ^ Eto W + !) - a (^o + + 2) 2 - (Z + h + 2) 2 (for details see [Hj). Note 
that the relation r y + 5 + e = a + P + l is satisfied. 

Equation (jBTf) with the condition 7 + 5 + e = a + j3 + 1 is called the Heun equation. 
The Heun equation is the standard canonical form of a Fuchsian equation with four 
singularities. By the aforementioned correspondence between x and w, the parameters 
(h, h, hi hi E, r) correspond essentially one-to-one to the parameters (a, j3, 7, 5, e, g, t) 
with the condition that 7 + 5 + e = a + /3 + 1 (see [TTH IT^l IT7| ) . 

Remark. The expression of the Heun equation in terms of elliptic functions as in 
(JMj) was essentially established by Darboux more than 100 years ago (see [T^fj). 

There are other possible choices for the function <3>(u>) of (j6f>)) and the relation (jfifij) 
between x and w. For example, Smirnov ^3] chose the correspondence defined by 



t = «aza = w -h/2 [w 

e 2 -ei' V / 

The transformation w 



(w - iy h/2 {w - t)- lz / 2 , and 



w 



(x)—ei 
ea-ei 



ei— e3 
p(x)-e 3 



produces the following correspondence: 



(68) 



a; 





1 
2 


T+l 


T 

2 


u> 





1 


-f- 

O, 


OO 



Other mathematical objects and concepts discussed in this paper are transformed 
under (JBTIj) as follows: 

• In sections ED and El we considered the finite-gap property for the Inozemtsev model 
for the case h,h,h,h G Z. The condition hih,h,h £ Z is transformed into the 
condition 7, 5, e, a — [3 <G | + Z. 

• Some statements in this paper have assumptions such as U = 0. If the assumption 
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Figure 1. 



li = is changed to Zj = — 1, all the applicable statements in this paper remain valid. 
The condition ! = or —1 (resp. l\ = or —1, l 2 = or —1, l 3 = or —1) is 
transformed under into the condition 7 = | or | (resp. <5 = |or|,e = |or|, 
a-/3 = ±|). 

• In section 12721 we introduced the doubly periodic functions up to signs. The condition 
of doubly periodicity up to signs in x is transformed into algebraicity in w. 

• By translating the functions (e.g., A(x, E), E(x, E)) and the formulas (e.g., (jSJ), (JBJ 
and ([211) ) m this paper or in [TB] into a form appropriate for the Heun equation (|67p. 
we recover the corresponding formulas of Smirnov |14j . 

• Let C be the cycle enclosing two points and 1 in w. The path from e to e + 1 (|e| 
sufficiently small) in x corresponds to the cycle C in w (see Figure 1). 

Assume 7, 8, e, a — (5 G \ +Z and 7 + 5 + e = a + /3 + 1. The coefficients of the mon- 
odromy on the cycle C are written as the hyperelliptic integral in q (the parameter 
defined after equation (jHTjl ). Concretely let Ai(w,g) and A 2 (u>,g) be the functions 

defined by the relations A, (^, -± + c )) £ (^) = A(s, £) and 

A 2 (#F^'-il(^ + c o))^(^) 7 A(-^^)- Then the functions A^w, g) 
and A 2 (w,g) satisfy equation (|57jl . Let Q(q) and Qi(g) be the monic polynomi- 
als defined by Q(E) = (4(e 3 - e 2 )) 2 * +1 Q + c )) and Q 1 (E) = (4(e 3 - 

e 2 )) 9 Qi {7^, + c o))> fi x Q* an d assume that Q(g*) 7^ 0. Finally let Af(w,g„) 

and A 2 (w, g*) be the functions obtained by analytic continuation in w along the cycle 
C from the functions Ai(w, g#) and A 2 (u>, g*), and define the monodromy matrix M c 
by 

(69) (A?(w,g,),Ajj(w,gv)) = (Ai(w, g*), A 2 (tw, g*))M c . 



Let g be a solution to the equation Q(qo) = 0. Then we have Ai(w, go) = A 2 (w, g ) 
and A^(w,g ) = ±A 1 (w,q ) (c.f. Corollary 13. 8|) . From formula (J25j) M c is a diagonal 




HEUN EQUATION III 



23 



matrix and the eigenvalues of M c are given by the following functions: 
(70) ± (-1)^ exp I V^T^J r 4==dq) , 

V Jqo yJ-Q{q) J 

± (-l) h+h exp f - y/te^Ttj f 4= d( H > 

V V-Q(?) / 

where the sign ± is chosen to be the same as the one A^w, go) = ±Ai(w, q ). Note 
that the factor (— 1)'°+*! comes from the branch of the function $(w). 

For the case Q(q») = it is necessary to select the other basis for the solutions to 
(J67j) which entails other calculations. 



Appendix A. Elliptic functions 

This appendix presents the definitions of and formulas for elliptic functions. 

Let u>i and 003 be complex numbers such that the value 003/001 is an element of the 
upper half plane. In this paper we set Ui — 1/2 and oo 3 = r/2. 

The Weierstrass p-function, the Weierstrass sigma-f unction and the Weierstrass 
zeta-function with periods {2uoi,2uo 3 ) are defined as follows: 



^ ^ ^ ^ z 2 ( (z — 2mwi 



(m,n)GZxZ\{(0,0)} 

2 



2nco>3) 2 (2mco>i + 2nus) 2 / 



a{z) =z J] ( 1 

(m,n)eZxZ\{(0,0)} 



2wjj0\ + 2nco>3 
2 



z 2 
■ exp I h 



2miO\ + 2nu>3 2(2mwi + 2na> 3 ) 2 / ' 



„ , a'[z) 
C{z) 



a{z) 

Setting 002 = —ooi — 003 and 

(72) e; = p(^), 77^ = C(oOi) (i = 1, 2, 3) 
yields the relations 

(73) e 1 + e 2 + e 3 =r] 1 +r]2 + r] 3 = 0, 

= -£(x), (p'(z)) 2 = A(p(z) - ei )(p(*) - ea)(p(a() - e 3 ), 

p(z + 2u j ) = fp(z), C(z + 2u j ) = C(z) + 2r lj (j = 1,2,3), 

/ \ (cj — ej/)(ej — ej») 
^ + ^ = e, + ^ ^ ^ z = 1, 2, 3 , 

where G {1,2,3} with i! < i", i 7^ i', and i 7^ z". The constants #2 and #3 are 
defined by 

(74) g 2 = -4(e x e 2 + e 2 e 3 + e 3 ei), 5-3 = 4eie 2 e 3 . 
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The co-sigma functions <Ti(z) (i = 1,2,3) and co-p functions pi{z) (i = 1,2,3) are 
defined by 
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